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Abstract 

Since the birth of the quark model, the diquark which is composed of two quarks 
has been considered as a substantial structure of color anti-triplet. This is not only 
a mathematical simplification for dealing with baryons, but also provides a physical 
picture where the diquark would behave as a whole object. It is natural to ask 
whether such a structure is sufficiently stable against external disturbance. The 
mass spectra of the ground states of the scalar and axial-vector diquarks which are 
composed of two-light (L-L), one-light-one-heavy (H-L) and two-heavy quarks (H- 
H) respectively have been calculated in terms of the QCD sum rules. We suggest a 
criterion as the quantitative standard for the stability of the diquark. It is the gap 
between the masses of the diquark and ^/sq where sq is the threshold of the excited 
states and continuity, namely the larger the gap is, the more stable the diquark 
would be. In this work, we calculate the masses of the type H-H to complete the 
series of the spectra of the ground state diquarks. However, as the criterion being 
taken, we find that all the gaps for the various diquaks are within a small range, 
especially the gap for the diquark with two heavy quarks which is believed to be 
a stable structure, is slightly smaller than that for other two types of diquarks, 
therefore we conclude that because of the large theoretical uncertainty, we cannot 
use the numerical results obtained with the QCD sum rules to assess the stability 
of diquarks, but need to invoke other theoretical framework. 
PACS numbers: 11.55.Hx, 12.38.Lg 



1 Introduction 

Right after the birth of the quark model, the diquark model was proposed: two quarks 
constitute a color-anti-triplet which may be a tightly bound state. In Gell-Mann's pioneer 
paper on the quark model, he discussed the possibility of existence of free diquarks jl]. 
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With the diquark picture, numerous authors have studied the processes where baryons 
are involved [Zl [3j SI El E] and their conclusions support the existence of diquarks. Even 
in the meson sector, some newly observed resonances are considered to possess the exotic 
structures. One possibility is that the mesons are tetraquarks composed of a diquark and 
an anti-diquark, [7J [8j |9j [TOj E]. In fact, it is still in dispute that the diquark is a real 
spatially bound state as a pseudo-particle or just a loosely bound state. Recently, the 
authors of Refs. [T2l IT3] treat the diquark as an explicit particle which is the essential 
ingredient inside hadrons (baryons or exotic mesons). If these diquarks indeed exist as 
stable particles, they should have certain and definite masses and quantum numbers. Just 
as we discuss the regular hadrons, their spectra not only possess the real part, but also the 
imaginary part which corresponds to the stability of the diquark. Namely, the lifetime 
of the diquark might be finite. The main scenario to determine the diquark lifetime 
is that with an external disturbance, the diquark might dissolve into two quarks. By 
our intuition, the diquark composed of two light quarks might be easier to dissolve by 
absorbing gluons. Generally, it is believed that the heavy diquarks which are composed 
of two heavy quarks are more stable against external disturbance. 

In principle, a baryon which is composed of three valence quarks is described by the 
Faddeev equation group composed of three coupled differential equations jTH [15]. For the 
equations, the three valence quarks are of the same weights. It is possible that two of the 
three quarks would accidentally constitute a bound state, say, by quantum fluctuation. 
The diquark can be treated as a sub-system and behaves as an independent particle. 
In this picture the baryon possesses a diquark-quark structure and thus the three-body 
system turns into a two-body one. Correspondingly, the three Faddeev equations reduce 
to a single equation (no mater relativistic or non-relativistic). Therefore the problem 
is greatly simplified. It is noted on other side, that the sub-system is not exactly a 
fundamental one, but possesses an inner structure. When it interacts with gluons, a form 
factor which manifests the inner structure should be introduced. Now, we are confronting 
a problem: which two valence quarks of the three would tend to combine into a bound 
state. In fact, any two quarks have a chance to combine via strong interactions, but the 
rest valence quark would interact with the individual quarks in the diquark and tend to 
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tear it apart. Thus the key point is if such sub-system is sufficiently stable against the 
disturbance. When we deal with the baryons which are composed of three light quarks, 
one-heavy-two-light quarks and two-heavy-one-light quarks, we notice an obvious unequal 
structures. We need to determine which type of diquark i.e. the diquarks with light-light 
(L-L), heavy-light (H-L) or heavy-heavy (H-H) structures, is more stable. Thus we can 
more confidently reduce the three Faddeev equations to one with the expected diquark 
subsystem. For serving this aim, it is significant to investigate the stability of diquarks. 

First of all, it is important and interesting to investigate the spectra of the diquarks 
with various quark contents. According to the masses of different flavors, one can cate- 
gorize the diquarks into three types: light-light diquark (L-L), light-heavy diquark (H-L) 
and heavy-heavy diquark (H-H), where light quarks are u, d, s and heavy quarks are c, b. 

The QCD is by all means a successful theory for the strong interaction, but the non- 
perturbative QCD which dominates the low energy phenomena is still not fully understood 
yet. Among the theoretical methods for treating the non-perturbative QCD effects, the 
QCD Sum Rules[16j is believed to be a powerful means for evaluating the hadronic spectra 
and other properties of hadrons. In Refs.|10[ITT]. the mass spectra of the scalar light-light 
diquark states (L-L) were studied with the QCD sum rules. Recently, with the same 
method, Wang studied the light-heavy diquark states (H-L)[18j. 

By the common sense, the diquark composed of two heavy quarks (H-H) may be a 
kinematically favorable sub-structure in baryon. For studying the stability of the three 
different types (L-L, H-L and H-H) diquarks, it would be crucial to discuss their properties 
in a unique framework. In this work, we are going to carry out the job in terms of the 
QCD sum rules. Then we will discuss the feasibility of such a scheme by scanning the 
numerical results obtained in this method. 

In this paper, we calculate the masses of heavy-heavy diquark states (H-H) with the 
QCD sum rules. Then combining the results presented in the relevant works [TTJ HH], we 
discuss the stability of the diquarks altogether. 

Obviously it is crucial to set a reasonable and practical criterion of the stability of the 
sub-structures — diquarks. 
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In the scenario of the QCD sum rules, numerically there exists a threshold s corre- 
sponding to a starting point beyond which higher excited and continuous states reside. 
This cutoff provides a natural criterion which we may use to study the stability of the 
diquark. Namely, according to our general knowledge on quantum mechanics, the con- 
tinuous spectra would correspond to the dissolved state where the constituents of the 
supposed-bound state would be set free. Thus we choose the gap between the ground 
state and the corresponding threshold t/so as the criterion of stability of the diquarks. 

However, as well known that there is a 20% theoretical uncertainty in all the com- 
putations via the QCD sum rules, therefore, even though such criterion may be indeed 
reasonable in principle, it is still doubtful if the results obtained in terms of the QCD 
sum rules can practically apply to reflect the diquark stability. The goal of this work is to 
testify the reasonability of applying the supposed criterion within the framework of the 
QCD sum rules. 

In the last section, we will come back to discuss the feasiability based on the numerical 
results we obtained in term of the QCD sum rules. 

The paper is organized as follows. After the introduction, in Sec. II we derive the 
correlation function of the suitable currents with proper quantum numbers in terms of 
the QCD sum rules. In Sec. Ill, our numerical results and relevant figures are presented. 
Section IV is devoted to a summary and concluding remarks. The tedious analytical 
results are collected in the appendices. 

2 Formalism 

For studying the scalar and axial-vector H-H diquark with the QCD sum rules, we 
write down the local color anti-triplet diquark currents: 



where k are the color indexes, Q — b,c, and C is the charge conjugation operator. 




e^ fc Qj(x)C 7 5 Q fc (x) 
e^ k Qj(x)C^Q k (x) 



(1) 
(2) 
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In order to perform the QCD sum rules, we define the two-point correlation functions 
11(g) (for scalar diquark) and 11^ (g) (for axial- vector diquark) as follows: 

11(g) = i J d l xe i * x (0\T{j i (x)J i \0)}\0) , (3) 

n^Cff) = z|d 4 xe^(o|r{j;(x)j:t( )}|o). (4) 

On the hadron side, after separating out the ground state contribution from the pole 
terms, the correlation function is expressed as a dispersion integral over the physical 
regime, 

n(?) = T7^ + ~ Hds^l, (5) 
M 2 -q 2 7T J s o s-q 2 



(6) 



where M t with subscript t being S or A for the scalar or axial-vector respectively, is 
the mass of the ground state diquark, p$(s) is the spectral density and represents the 
contribution from the higher excited states and the continuum, s° is the threshold for the 
excited states and continuum, the pole residues A f correspond to the diquark coupling 
strength, is defined through [IB]: 

(0|J*(0)|^(g)) = X s S ij , (7) 
(0|Jj(0)|^(g)> = Xas^, (8) 

with e M being the polarization vector of the axial-vector diquark. 

On the quark side, the operator product expansion (OPE) is applied to derive the 
correlation functions. Firstly, we can write down the "full" propagator Sq (x) of a massive 
quark, where the vacuum condensates are clearly displayed [19]. 



WW = / ***** 



8a 9s{t k ) ij Gf a af} tf + m Q ) + (jt + m Q )a aP 



QK ' ' (2tt) 4 J F \$-m Q 4 (p 2 -m 2 Q ) 2 

it 2 a s GG p 2 + m Q p : . 

+ T ( ~ ) ^ mq (^-m^ + -^ (9) 



5 



where ( asGg ) = { asGa P G — ^ then contracting the quark fields in the correlation functions, 
we gain the results: 

n(g) = -ie i] h ij ' k ' J d 4 xe iq - x Tr!^ 5 S^\x)^CS^' T (x)C^ , (10) 

IV(g) = ie^ k e^' k ' J d 4 xe iq x Tr^S^\x) lu CS kk ' T (x)cY (11) 

Then substituting the full c and 6 quark propagators into above correlation functions 
and integrating over the variable k, we obtain the correlation functions at the level of 
quark-gluon degrees of freedom. Simply, the correlation function n t (g 2 ) (t=S or A) is 
written as: 

U t (q) = Ur t (q) + Iir dA (q), (12) 

where the superscripts "pert", and"cond" refer to the contribution from the perturbative 
QCD, and gluon condensates, respectively. In this work, we only keep the two-gluon 
condensate in consideration for the heavy quark condensates are zero as suggested in 
literature [19]. 



Due to the quark-hadron duality, we differentiate Eq. tTHl) with respect to y^-, then 
eliminate the pole residues At, and obtain the resultant sum rule for the mass spectra of 
the H-H diquark states: 



M *=\l-%> ( 13 ) 



with 



R° t = - p dspr\s)e- s / M B + B[Ur dA (q 2 )}, (14) 

R\ = tt~~~ — o-Rt ■ (15) 

Here, Mb is the Borel parameter and Sq is the threshold cutoff introduced to remove the 
contribution of the higher excited and continuum states [201 l21j . 

For the H-H diquark states, the detailed expressions of are collected in the appendix. 
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3 Numerical Analysis 



3.1 The masses of the ground diquark states with only heavy 
flavors. 

The numerical parameters used as inputs in this work are taken as [T6l [T8l [I~9] 122] 

(qq) = -(0.24 ± O.OlGeV) 3 , (ss) = (0.08 ± 0.02) (qq) , 

(qg s (rGq) = ml(qq) , (sg s aGs) = m 2 (ss) , 

m 2 = (0.8 ± 0.2)GeV 2 , (-G 2 ) = 0.012GeV 4 , (16) 

m u ~m d = 0.005GeV , m s = (0.14 ± 0.01)GeV , 

m c = (1.35 ± 0.10)GeV , m b = (4.7 ± 0.1)GeV . 
where the energy scale is \i = 1 GeV. 

It is crucial to determine the proper threshold s° and Borel parameter M\ for obtaining 
physical spectra. For justifying if the choice is suitable, there are two criteria. First, 
the perturbative contribution should be larger than the contributions from all kinds of 
condensates, and another is that the pole contribution should be larger than the continuum 
contribution [16j [19]. In our work, the error bars are estimated by varying the Borel 
parameters, s° and including the uncertainties of the input parameters as well. 





mass (GeV) 


M2(GeV 2 ) 


v^GeV) 


pole 




cc(l+) 


2.99 ±0.10 


1.2 - 2.5 


3.3 ±0.1 


(56 - 88)% 


(4 - 8)% 


bc(0+) 


6.30 ±0.09 


2.2 - 5.0 


6.6 ±0.1 


(85 - 99)% 


(5 - 14)% 


6c(l+) 


6.36 ±0.08 


3.0-6.0 


6.7±0.1 


(53 - 85)% 


(3 - 5)% 


bb(l + ) 


9.76 ±0.08 


6.0 - 14.0 


10.1 ±0.1 


(41 - 79)% 


(0.04 - 0.23)% 



Table 1: For the H-H diquark states, we show the masses, the preferred Borel parameters 
Mg, the threshold parameters s , the contribution from the pole term to the spectral 
density and the contribution from {^-G 2 ). 

The spectra of the L-L and H-L diquarks have been studied in the previous literature [T7] 
ITS] . Our input parameters are the same as that used in Ref.[TB]. Numerically we have 
rechecked the results of the previous works for spectra of L-L and H-L diquarks and find 
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Figure 1: Dependence of M^ c , Mjf, M| c and Mjf on the Borel parameter M|. We 
deliberately put two vertical lines denoting the chosen Borel window. 



our results are consistent with them, thus we simply present the relevant results for the 
L-L and H-L diquarks in the following table. The analytical formulations for the H-H di- 
quarks are shown in the appendix and that for L-L and H-L can be found in the relevant 
references which we list in the bibliography of the work. 



3.2 Discussion on the stability of the ground diquark states. 

In Table[ll we present the masses and the threshold ^/sq of the H-H scalar and axial- 
vector diquark ground states. Then, altogether with that given in literature, we collect 
the results for all these three types of diquark in the following table. Here we define the 
energy gap as AE = ^/so — M d , where M d is the mass of the corresponding diquark state. 

As aforementioned, we use energy gap AE to embody the stability of diquarks. Namely, 
according to our general knowledge on quantum mechanics, the continuous spectra would 
correspond to a dissolved state where the constituents of the supposed-bound state are set 
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mass(GeV) 


M B (GeV 2 ) 


Vs°(GeV) 


AE(GeV) 


sq(0 + ) 


0.55 ±0.03 


1.0 


0.95 


0.40 


qq(l+) 


0.34 ±0.04 


4.0 


0.85 


0.51 


sq(l+) 


0.42 ±0.03 


4.0 


0.95 


0.53 


ss(l + ) 


0.50 ±0.05 


4.0 


1.05 


0.55 



Table 2: The diquark masses, the preferred Borel parameters M#, the preferred threshold 
parameters s° obtained by choosing reasonable plateaus, and the energy gap AE for each 
possible L-L diquark state. Here the revelent results are derived from the corresponding 
formulae in Ref.[T7] with our input parameters. 





mass(GeV) 


M^GeV 2 ) 


v^°(GeV) 


A£(GeV) 


cg(0+) 


1.77 


1.50 


2.19 


0.42 


cg(l+) 


1.76 


1.60 


2.19 


0.43 


cs(0+) 


1.84 


1.55 


2.24 


0.40 


cs(l+) 


1.84 


1.65 


2.24 


0.40 


bq(0 + ) 


5.14 


3.85 


5.48 


0.34 


bq(l+) 


5.13 


3.95 


5.48 


0.35 


bs(0+) 


5.20 


3.95 


5.57 


0.37 


bs(l+) 


5.22 


4.10 


5.57 


0.35 



Table 3: The diquark masses, the preferred Borel parameters Mb, the preferred threshold 
parameters s° obtained by choosing reasonable plateaus for the H-L diquarks. These 
revelent results have already been obtained in Ref.|18j. 

free. Thus we choose the gap between the mass of the ground state and the corresponding 
threshold y/so as the criterion of stability of the diquark. In Tables([2]|4]), we show the 
energy gap AE for the three types of diquarks. Concretely, the energy gap is between 0.40 
GeV and 0.55 GeV for L-L diquarks, between 0.34 GeV and 0.43 GeV for H-L diquarks, 
and between 0.30 GeV and 0.34 GeV for H-H diqaurks. One can observe that the sequence 
of the energy gaps for these three types of diquarks is :AEl-l > AEh-l > AEh-h- If 
larger energy gap implies more stable structure, the L-L diquarks should be the stablest 
one and the H-H diquarks is the most instable structure, however, this definitely contra- 
dicts to our intuition. Let us discuss this issue in next section. 
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mass(GeV) 


MJ>(GeV 2 ) 


Vs°(GeV) 


AS(GeV) 


£>c(0+) 


6.30 


3.6 


6.6 


0.30 


cc(l+) 


2.99 


1.9 


3.3 


0.31 


6c(l+) 


6.36 


4.5 


6.7 


0.34 


66(1+) 


9.76 


10.0 


10.1 


0.34 



Table 4: The H-H diquark masses, the preferred Borel parameters Mg, the preferred 
threshold parameters s° obtained by choosing reasonable plateaus. 

4 Summary and Conclusions 

In this work, we try to study the stability of scalar and axial vector diquarks of three 
types L-L, H-L and H-H in terms of the QCD sum rules. For the purpose, we first calculate 
the mass spectra of the scalar (0 + ) and axial- vector (1 + ) H-H diquarks. Together with the 
results given in the previous works about the L-L diquarks [T7] and the H-L diquarks |1 8], 
we can compare the corresponding AE which is supposed to be a reference criterion for 
the stability of the diquark sub-structure. 

However, the numerical results obtained in terms of the QCD sum rules determine the 
sequence AE L _ L > AE H _ L > AE H ^ H . Even though actually the values of AE L _ L , AE H _ L 
and AEl_l are not far apart, this sequence obviously conflicts to our common sense where 
the H-H diquark is the most stable structure. 

Because the computation with the QCD sum rules possess over 20% errors and the 
determined AE values do not differ much (less than 20%), we cannot decide that the H-H 
diquark is the most instable one. 

Thus, we would be tempted to conclude that due to the large errors brought up in 
the computations with the QCD sum rules, it is not suitable to determine the stability 
in terms of the QCD sum rules, and we have to turn to other theoretical approaches to 
investigate this important issue which would be the goal of our next work. 
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Appendix 

For the H-H diquark, our analytical expressions are shown as follows: 
3 



Ps{s) 



4ns 



(m 2 Ql + 2m Ql m Q2 + m 2 Q2 - s) yj (m 2 Qi - m\ 2 + s) 2 - Am 2 Qi s , (17) 



Pa{s) = ^ {m 2 Qi - Am Ql m Q2 + mJJ yj (m 2 Qi - m 2 Q2 + s) 2 - Am 2 Ql s , (18) 



GUM! 



G 2 (M 2 B ) 



2 2 
Ql i Q2 



(a s G 2 ) / dxe 
Jo 



(4-9x)(l-x) 
16tt 



+ 



32vrx 3 M| 



2(1 -x)(9x z 
1 

+ 



-13x + A)x 3 m 2 Qi + 18xm Ql m Q2 + 2{x - l)(9x - A)x A m 2 Q ^ -r 

x (1 - x)(-9x 5 + 22a; 4 - 17a; 3 + Ax 2 -7x + 7)xm 4 Ql + (6x 3 - Qx 2 + 3x 
-2)m^m Q2 + {l-x) (18a; 4 - 26a; 3 + 8a; 2 + 7) x 2 m 2 Q m 2 Q2 - 6x 3 m Ql m 3 Q2 



1 

32irx 3 M% 



3(l-x)(a:-l) 2 (a; 3 -a; 2 + l) 



+ (4 - 9a;) (1 - x)x 5 rriQ 2 

xm 6 Ql + 3 (a; 4 - 2a; 3 + x 2 + x - l) xm 5 Ql m Q2 + 3(x - l)(3x 4 - 6a; 3 + 
+2x - 2)x 3 rriQ i m 2 Q2 - 3 (2x 3 - 2a; 2 + l) x 2 m 3 Q m 3 Q2 + (1 - x) (9a; 3 - 9a; 2 



3x z 



+3)x 4 m 2 3i m 4 2 + 3a; 5 m Ql m^ 2 + 3(x - l)a; 7 mj 2 
Gf(Mi,m 0l ^m Q2 ) , 



(19a) 
(19b) 
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m Qi , m Q2 



G 3 (M 2 B ) = (a s G 2 } J dxe * M ^ { 



9 (2x 2 - 2x - 1) 



x 



64tt 64tt(x - l) 3 x 3 M 2 

3(3x 2 (2x 2 - 2x - 1) (x - lfm 2 Ql + 8a; 2 (x - l) 2 m Ql m Q2 

3(3(z - lfx 4 m 4 Q2 + 3(x - 1) 4 iV 4 0i + (x - l)x 2 (-6a: 3 + 12x 2 



x 



-7x+ l)m 2 Q2 m 2 Ql ] 



3(-x b m b Q2 + (x- l)x° 



64n{x - l) 3 x 3 M|, 

(3x - 2)m% 2 m 2 Ql + (x - l) 5 m 6 Qi + (1 - 3x)(x - l) 3 xm 2 Q2 m 4 Ql )] } .(20a) 



G\(Mj 



2 2 

m Qi , m Q 2 



{a s G 2 ) / dxe 



(x- l)(9x-4) 



-2x d 



k 8tt 1 (iv: (./- 1)V\\/ 2 

x(9x-4)(a;- l) 2 mg 2 + 2x 2 (9x - 4)(x - l) 3 m 2 Ql + 36x(x - l) 2 m Ql m Q 
1 

~167r(a;- l) 2 x 3 M| 
-(a; - l) 2 x (I8x 2 - 8x - 7) m 2 Q2 m 2 Qi + 12(x - l)x 2 m 3 Q2 m Ql - 2(x - 1) 



[x - l)x 3 {9x - A)m 4 Q2 + {x- l) 3 (9x 2 - 4x - 7) m 4 Ql 



x(3x + 2)m Q2 m 3 Qi 



16n(x - l) 2 x 3 M B 



- 3x 4 m 6 Q2 + 6x 6 m° Q2 m Ql 



3^5 



+(x - l)x 2 (9x - 3)m% 2 m 2 Qi + 3{x - l) 4 m 6 Qi + 6(x - l) 3 m Q2 m^ 
-3(x - l) 2 x(3x - 2)m 2 Q2 m 4 Ql - 6(x - l)x(2x - l)m 3 Q2 m 3 Ql ] } , (21a) 
G 2 A (M 2 B ) = G?(M 2 B ,m Ql ^m Q2 ), (21b) 
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G\{M%) 



/l x ~^ 1 x ( Q 

cfee "5 | - - 





x 3 ((x - l) 3 x 2 (6x 2 - 6x + 5) m 2 Qi - (x - l) 2 x 3 (6x 2 - 6x + 5) m| 2 ) 



(22a) 
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